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Abstract 

We calculate the two-loop QCD correction to the scalar quark contributions to the 
electroweak gauge boson self-energies at zero momentum-transfer in the supersymmetric 
extension of the Standard Model. We then derive the 0(a s ) correction to the contribution 
of the scalar top and bottom quark loops to the p parameter, which is the most sizable 
supersymmetric contribution to the electroweak mixing angle and the W^-boson mass. The 
two-loop corrections modify the one-loop contribution by up to 30%; the gluino decouples 
for large masses. Contrary to the SM case where the QCD corrections are negative and 
screen the one-loop value, the corresponding corrections in the supersymmetric case are 
in general positive, increasing the sensitivity in the search for scalar quarks through their 
virtual effects in high-precision electroweak observables. 
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Supersymmetric theories (SUSY) |J are the best motivated extensions of the Standard 
Model (SM) of the electroweak and strong interactions. They provide an elegant way to 
break the electroweak symmetry and to stabilize the huge hierarchy between the GUT and 
the Fermi scales, and allow for a consistent unification of the gauge coupling constants as 
well as a natural solution of the Dark Matter problem; for recent reviews see Ref. M. 

Supersymmetry predicts the existence of scalar partners to each SM fermion, and spin- 
1/2 partners to the gauge and Higgs bosons. So far, the direct search of SUSY particles 
at colliders has not been successful, and under some assumptions one can only set lower 
bounds of (9(100) GeV on their masses [Q]. The search can be extended to slightly larger 
values at LEP2 and the upgraded Tevatron; higher energy hadron or e + e~ colliders will be 
required to sweep the entire range of the SUSY particle masses up to the TeV scale. 

An alternative way to probe SUSY is to search for the virtual effects of the additional 
particles. Indeed, now that the top-quark mass — the measured value of which being 
in remarkable agreement with the predicted one — is known 0], one can use the high- 
precision electroweak data to search for the quantum effects of the SUSY particles: sfermions, 
charginos/neutralinos and gluinos. 

In the Minimal Supersymmetric Standard Model (MSSM) it is well known that, besides 
the rare decay b — > sj [Q, there are two possibilities for the virtual effects of SUSY particles 
to be large enough to be detected in present high-precision experiments. The first possibility 
is that charginos and scalar top quarks are light enough to affect the decay width of the 
Z boson into 6-quarks |J; however, for masses beyond the LEP2 or Tevatron reach, these 
effects become too small to be observable ||. 

The second possibility is the contribution of the scalar top and bottom quark loops 
to the electroweak gauge-boson self-energies 0: if there is a large splitting between the 
masses of these particles, the contribution will grow with the mass of the heaviest scalar 
quark and can be sizable. This is similar to the SM case, where the top/bottom weak 
isodoublet generates a quantum correction that grows as mf. This contribution enters the 
electroweak observables via the p parameter ||, which measures the relative strength of 
the neutral to charged current processes at zero momentum-transfer. It is mainly from this 
contribution that the top-quark mass has been successfully predicted from the measurement 
of the effective electroweak mixing angle at the Z-boson resonance and the P^-boson mass 
at hadron colliders ||, a triumph for the electroweak theory. 

In order to treat the SUSY loop contributions to the electroweak observables at the 
same level of accuracy as the standard contribution, higher order corrections should be 
incorporated. In particular the QCD corrections, which because of the large value of the 
strong coupling constant can be rather important, must be known. It is the purpose of this 
report to provide the two-loop QCD corrections to the scalar quark contributions to the 
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electroweak precision observables. As a first step, we will consider here only the contributions 
to the p parameter; more detailed results will be given elsewhere 0. 

The p parameter, in terms of the transverse parts of the W— and Z-boson self-energies 
at zero momentum-transfer, is given by 

i a _ n zz (o) rwg 

9 1 - Ap ' * P M 2 Myf [L) 

In the SM, the contribution of a fermion isodoublet (u, d) to Ap reads at one-loop order 

A,r^F„Km 2 d ) , (2) 

with the color factor N c and the function Fq given by 

Fo(x,y) = x + y- — ^-log- . (3) 
x-y y 

The function F vanishes if the u- and rf-type quarks are degenerate in mass: F (m 2 , m 2 ) = 
0; in the limit of large quark mass splitting it becomes proportional to the heavy quark mass 
squared: F (m 2 ,0) = m 2 . Therefore, in the SM the only relevant contribution is due to 
the top/bottom weak isodoublet. Because m t 3> m b , one obtains Ap^ M = 3G pm 2 / (8\/2ir 2 ) , 
a large contribution which allowed for the prediction of m t . However, in order that the 
predicted value agrees with the experimental one, QCD corrections have to be included. 



These two-loop corrections have been calculated ten years ago, leading to a result |T0 
Apf M = -Apf* ■ + tt 2 /3). For the value a s ~ 0.12, the QCD correction fly] decreases 
the one-loop result by approximately 10% and shifts m t upwards by an amount of ~ 10 GeV. 

In SUSY theories, the scalar partners of each SM quark will induce additional contribu- 
tions. The current eigenstates, qi and (Jr, mix to give the mass eigenstates. The mixing 
angle is proportional to the quark mass and therefore is important only in the case of the 



third generation scalar quarks ||12|| . In particular, due to the large value of mt, the mixing 
angle 0j between t L and in can be very large and lead to a scalar top quark t\ much lighter 
than the t-quark and all the scalar partners of the light quarks [|1^]. The mixing in the 
6-quark sector can be sizable only in a small area of the SUSY parameter space. 

The contribution of a scalar quark doublet (u, d) to the transverse parts of the W/ Z- 
boson self-energies at zero momentum-transfer [Fig. 1] can be written as |7j 



nw(o) = -T* E («>|) 2 *b« )W 

V^TT i,j=l,2 



■2 



V (7—u.fi 



q — u,d 



where the factors a\ are given in terms of the scalar quark mixing angle 9g as a\ = cosOg 
and a\ = sin 9g. 
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As can be seen from Fq in eq. (3), the contribution of a scalar quark doublet vanishes if all 
masses are degenerate. This means that in most SUSY scenarios, where the scalar partners 
of the light quarks are almost mass degenerate, only the third generation will contribute. 
Neglecting the mixing in the b sector, Ap is given at one-loop order by the simple expression 

+ cos 2 e~ t F (m| , mfj + sin 2 9~ t F (m| , mf J] . (5) 

In a large area of the parameter space, the scalar top mixing angle is either very small 6 t ~ 
or maximal, 9 t ~ — 7r/4. The contribution A/)q USY is shown in Fig. 2 as a function of the 
common scalar mass = mz = m~ h for these two scenarios. The contribution can be at 
the level of a few per mille and therefore within the range of the experimental observability. 
Relaxing the assumption of a common scalar quark mass, the corrections can become even 
larger 0. 

At 0(aa s ), the two-loop Feynman diagrams contributing to the p parameter in SUSY 
[Fig. 3] consist of two sets which, at vanishing external momentum and after the inclusion 
of the counterterms, are separately ultraviolet finite and gauge- invariant. The first one has 
diagrams involving only gluon exchange, Fig. 3a; in this case the calculation is similar to 
the SM, although technically more complicated due to the larger number of diagrams and 
the presence of q mixing. The diagrams involving the quartic scalar-quark interaction in 
Fig. 3a will either contribute only to the longitudinal component of the self-energies or can 
be absorbed into the q mass and mixing angle renormalization as will be discussed later. The 
second set consists of diagrams involving scalar quarks, gluinos as well as quarks, Fig. 3b; 
in this case the calculation becomes very complicated due to the even larger number of 
diagrams and to the presence of up to 5 particles with different masses in the loops. 

We have calculated the two-loop contribution of a complete quark/squark generation 
to the vacuum polarization functions of the electroweak gauge bosons at zero momentum- 
transfer, taking into account general mixing between scalar quarks and allowing for all 
particles to have different masses. In the following, we summarize the main features of the 
calculation [ IB ]. 



Our results have been derived by two independent calculations using different methods. 
In one method, the unrenormalized self-energies together with the mass and mixing angle 
counterterms were calculated with the help of the program Process Diagram [JH |, while in 



the other case the packages FeynArts |T5| [in which the relevant part of the MSSM has been 



implemented] and Two Calc [|16| were used to generate and evaluate the full set of Feynman 
diagrams and counterterms. The two independent calculations allowed for thorough checks 
of the final results. 

The two-loop Feynman diagrams of Fig. 3 have to be supplemented by the corresponding 
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counterterm insertions into the one-loop diagrams. By virtue of the Ward identity, the vertex 
and wave-function renormalization constants cancel each other. The mass renormalization 
has been performed in the on-shell scheme, where the mass is defined as the pole of the 
propagator. The mixing angle renormalization is performed in such a way that all transitions 
from qi <-> q~j which do not depend on the loop-momenta in the two-loop diagrams are 



canceled; this renormalization condition is equivalent to the one proposed in Ref. [17] for 
scalar quark decays. With this choice of the mass and mixing angle renormalization, the pure 
scalar quark diagrams in Fig. 3a that contribute to the transverse parts of the gauge-boson 
self-energies are canceled. 

In order to discuss our results, let us first concentrate on the contribution of the gluonic 
corrections, Fig. 3a, and the corresponding counterterms. At the two-loop level, the results 
for the electroweak gauge boson self-energies at zero momentum-transfer have very simple 
analytical expressions. In the case of an isodoublet (u, d) where general mixing is allowed, 
the structure is similar to eq. (4) with the a\ as given previously: 

G F M^a s ^ / q j\2 / 
IW(0) = y= - ^ {a.aA F 1 (m Ui ,m s 



4\/27r J 

nzz(0) = "^tI? E «4) 2 * R <) • (6) 



q — u,d 
i,i'=l,2 



The two-loop function Fi(x,y) is given in terms of dilogarithms by 



F 1 (x,y)=x + y- 2—^— log - 
x-y y 



x x 
2 + - log - 

y y 



+ £±*£]og"5-2(«-»)U,(l-=Y (7) 
(x-y) 2 y V Vj 

This function is symmetric in the interchange of x and y. As in the case of the one-loop 
function F , it vanishes for degenerate masses, Fi(x,x) = 0, while in the case of large mass 
splitting it increases with the heavy scalar quark mass squared: Fi(x, 0) = x(l + vr 2 /3). 

From the previous expressions, the contribution of the (t, b) doublet to the p parameter, 
including the two-loop gluon exchange and pure scalar quark diagrams are obtained straight- 
forwardly. In the case where the b mixing is neglected, the SUSY two-loop contribution is 
given by an expression similar to eq. (5): 



A ^ USY = l-si^cos^ K 



Gf&s 



Ti" 



m 



t2 



AV2' 

+ cos 2 9iF x (m| , mf J + sin 2 e~ t F x (m| , m| J] . (8) 

The two-loop gluonic SUSY contribution to Ap is shown in Fig. 4 as a function of the 
common scalar mass m,q, for the two scenarios discussed previously: B~ t = and 6^ ~ — 7r/4. 
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As can be seen, the two-loop contribution is of the order of 10 to 15% of the one-loop result. 
Contrary to the SM case [and to many QCD corrections to electroweak processes in the SM, 
see Ref. for a review] where the two-loop correction screens the one-loop contribution, 
ApSUSY same s jg n as ApQ USY . For instance, in the case of degenerate t quarks with 

masses m t - ^> m^, the result is the same as the QCD correction to the (t, b) contribution in 
the SM, but with opposite sign. The gluonic correction to the contribution of scalar quarks 
to the p parameter will therefore enhance the sensitivity in the search of the virtual effects 
of scalar quarks in high-precision electroweak measurements. 

The analytical expressions of the contribution of the two-loop diagrams with gluino 
exchange, Fig. 3b, to the electroweak gauge boson self-energies are very complicated even 
at zero momentum-transfer. Besides the fact that the scalar quark mixing leads to a large 
number of contributing diagrams, this is mainly due to the presence of up to five particles 
with different masses in the loops. The lengthy expressions will be given elsewhere 0. It 
turned out that in general the gluino exchange diagrams give smaller contributions compared 
to gluon exchange. Only for gluino and scalar quark masses close to the experimental lower 
bounds they compete with the gluon exchange contributions. In this case, the gluon and 
gluino contributions add up to ~ 30% of the one-loop value for maximal mixing [Fig. 5]. 
For larger values of nig, the contribution decreases rapidly since the gluinos decouple for 
high masses. 

Finally, let us note that for the diagrams in Fig. 3a analytical expressions for arbitrary 
momentum-transfer can be obtained as will be discussed in Ref. 0. With the present 
computational knowledge of two-loop radiative corrections, analytical exact results for the 
diagrams involving gluino exchange, Fig. 3b, cannot be obtained for arbitrary q 2 \ either 
approximations like heavy mass expansions or numerical methods have to be applied. 

In summary, we have calculated the two-loop 0(a s ) correction to the scalar quark contri- 
butions to the weak gauge boson self-energies at zero momentum-transfer in SUSY theories, 
and derived the QCD correction to the p parameter. The gluonic corrections are of (9(10%): 
they are positive and increase the sensitivity in the search for scalar quarks through their 
virtual effects in high-precision electroweak observables. The gluino contributions are in gen- 
eral smaller except for relatively light gluinos and scalar quarks; the contribution vanishes 
for large gluino masses. The phenomenological implications of our results will be discussed 
in a forthcoming paper. 
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FIGURES 
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FIG. 1. Feynman diagrams for the contribution of scalar quark loops to the gauge boson 
self-energies at one-loop. 
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FIG. 2. One-loop contribution of the (t, b) doublet to Ap as a function of the common mass 
rriq, for 9f = and % ~ —tt/4 [with tan/3 = 1.6 and ttilr = and 200 GeV, respectively, where 
mLR is the off-diagonal term in the t mass matrix]. 
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FIG. 3. Typical Feynman diagrams for the contribution of scalar quarks and gluinos to the 
W/Z -boson self-energies at the two-loop level. 
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FIG. 4. Gluon exchange contribution to the p parameter at two-loop as a function of m,q for 
the scenarios of Fig. 2. 
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ne 5. Contribution of the gluino exchange diagrams to A/)^ USY for two values of nig in the 
scenarios of Fig. 2. 
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